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1 INTRODUCTION 



Observ ations by the EXOSAT and later by the Ginga satel- 



ABSTRACT 

Among the variability behaviour exhibited by neutron star systems are the so- 
called "horizontal branch oscillations" (HBO, with frequencies « 50 Hz), the "lower- 
frequency kHz quasi-periodic oscillation" (QPO) and the "upper-frequency kHz 
QPO" , with the latter two features being separated in frequency by an amount com- 
parable to, but varying slightly from, the suspected spin-frequency of the neutron 
star. Recently, Psaltis, Belloni, & van der Klis (1999) have suggested that there ex- 
ists a correlation between these three frequencies that, when certain identifications 
of variability features are made, even encompasses black hole sources. We consider 
this hypothesis by reanalyzing a set of GX 339—4 observations. The power spectral 
density (PSD) constructed from a composite of 7 separate, but very similar, obser- 
vations shows evidence for three broad peaks in the PSD. If the peak frequencies of 
these features are identified with "QPO" , then their frequencies approximately fit the 
correlations suggested by Psaltis, Belloni, & van der Klis (1999). We also reanalyze a 
Cyg X-l observation and show that the suggested QPO correlation may also hold, but 
that complications arise when the "QPOs" (which, in reality, are fairly broad features) 
are considered as a function of energy band. These fits suggest the existence of at least 
three separate, independent physical processes in the accretion flow, a hypothesis that 
is also supported by consideration of the Fourier frequency-dependent time lags and 
coherence function between variability in different energy bands. If these variability 
features have a common origin in neutron star and black hole systems, then "beat 
frequency models" of kHz QPO in neutron star systems are called into question. 

Key words: accretion — black hole physics — Stars: binaries — X-rays:Stars 



tures have been observed in lower-luminosity neutron star 
systems, spe cifically the so-called atoll sources (Btrohmayei 
et al. 1996). The atoll systems also have exhibited low- 



lites during 1980's and early 1990's, and more recently by the 
Rossi X-ray Timing Explorer (RXTE) have revealed a rich 
variety of X-ray variability behaviour in neutron star and 
black hole candidate (BHC) systems. Neutron sta r systems 
in general, and the s o-called Z-sources in specific (Hasinger 



& van der Klis 198£), have exhibited a set of relatively nar- 
row features in their X-ray variability power spectral den- 
sities (PSD). These features, referred to as quasi-periodic 
oscillations (QPO), phenomenologically appear to fall pre- 
dominantly into one of four classes. At low Fourier frequen- 
cies (see van der Klis 1995 for a review) there are the so- 
called normal branch oscillations (NBO) with / ~ 5-20 Hz 
and the horizontal branch oscillations (HBO) with / « 15- 
60 Hz. A second set of QPO, typically occurring in pairs at 
high frequencies / ~ 200 — 1200 Hz, are referred to as the 
"lower- frequency kHz QPO" an d the "upper-frequency kHz 
QPO" ( van der Klis et al. 1996| ) . Similar high frequency fea- 



frequency features that in some ways are similar to the HBO 
(see Hasinger & van der Klis 1989, Homan et al. 1998). 

A wide variety of features ranging from narrow (see, 
for example, Nowak, Wilms, & Dove 1999; hereafter NWD) 
to broad "noise components" (see, for example, the discus- 
sion of van der Klis 1994a, van der Klis 1994b) have been 
observed in BHC systems as well. Frequently, the PSDs asso- 
ciated with these systems are flat at low frequencies, show a 
low frequency break into a steeper / _1 -/ -2 spectrum, and 
often exhibit a low frequency QPO at frequencies slightly 
above the break (see Wijnands & van der Klis 1999, and 
references therein). It has been noted that the location of 
the br eak and the frequency of the QPO are often corre- 
lated (Wijnands & van der Klis 1999), both in neutron stars 
and BHC systems. 



In a recent work, Psaltis, Belloni, & van der Klis ( 199£ ) 
(hereafter PBK) have presented a further analogy between 
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the features observed in neutron star and in BHC systems. 
Specifically, they have shown that the HBO frequency is 
apparently correlated with the lower-frequency kHz QPO 
frequency (Fig. 1 of PBK) . With suitable identifications (in 
part relying upon the break freq uenc y- Q P O frequency cor- 
relation of Wijnands & van der Klis (1999), as discussed by 
PBK) to features observed in atoll and BHC systems, this 
correlation is seen to extend over three decades in Fourier 
frequency and encompass Z-sources, atoll sources, and BHC 
sources. If these apparent correlations exist because of an 
underlying common physical mechanism, then all these vari- 
ability phenomena are somehow intrinsic to the accretion 
flow, and do not explicitly rely upon the presence or ab- 
sence of either a hard surface or a magnetic field in these 
systems. 

A correlation is also observed between the lower- 
frequency kHz and upper- frequency kHz QPO. Some models 
(e.g., Miller, Lamb, & Psaltis 1998) associate this difference 
frequency with being nearly equal to the spin-frequency of 
the neutron star (in contrast to the theories of Psaltis & 
Norman , Stella & Vietri 1998, Merloni et al. 1999, and 
Stella et al. 1999 discussed further below). PBK, however, 
speculated about whether the observed lower-frequency kHz 
QPO/upper- frequency kHz QPO correlation would also ex- 
tend to low frequencies and BHC systems in a manner simi- 
lar to the putative HBO/lower-frequency kHz QPO correla- 
tion. In this work we reexamine a set of observations of the 
BHC GX 339—4 in order to search for evidence of a "third 
QPO" in this system. We present such evidence in In §|3j 
we also present evidence for multiple broad-peaked features 
in the PSD of Cyg X-l. We discuss in the significance 
of these fits as regards the hypothesized QPO correlations. 
Furthermore, in light of the observed time lags and degree 
of linear correlation (i.e., the coherence function) between 
the soft and hard X-ray variability, we discuss whether the 
individual "QPO fit-components" to the PSD actually rep- 
resent physically distinct processes in the accretion flow. We 
summarize our results in 80. 



2 POWER SPECTRAL DENSITIES FOR 
GX 339 4 

In a previous work (NWD), we presented timing analysis 
for eight separate RXTE observations of GX 339—4. The 
faintest observation, which in terms of 3-9 keV flux is a 
factor ~ 2.5 fainter than the next brightest observation, 
clearly shows larger amplitude and characteristically lower- 
frequency variability than the other observations. Further- 
more, the faintest observation shows time lags between its 
soft and hard X-ray variability (Miyamoto fc Kitamoto 19891; 



Miyamoto et al. 1992| ; [Vaughan fc Nowak 1997] ; |Nowak et al 



1999a) that are significantly shorter than those exhibited by 



the other observations. The remaining observations, which 
span a factor of ~ 2 in terms of 3-9 keV flux, have relatively 
similar timing properties. Specifically, the shapes and am- 
plitudes of their PSDs are all comparable, and each shows a 
narrow PSD peak between 0.26-0.34 Hz. In addition, their 
Fourier frequency-dependent time lags span less than a fac- 
tor of two at any given Fourier frequency. 

In the frequency range of 0.1-30 Hz, NWD were rela- 
tively successful in fitting the PSDs to a functional form 
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Figure 1. Fourier frequency times composite power spectral den- 
sity (PSD) of GX 339-4. Thick solid line is the best fit model, 
whereas the dashed lines are the separate components of the 
model. The lower panel shows the residuals from the fit. 



that consisted of a power law approximately tx / , 
two additional Lorentzian components of the form 



P(f) 



R 2 Qfo 



fS + Q 2 (f-fo) 2 



plus 



(1) 



where /o is the resonant frequency of the Lorentzian, Q 
is the quality factor (~ /o/A/, where A/ is the full- 
width-half-maximum of the Lorentzian), and R is the fit 
amplitude (root mean square variability, rms = -R[l/2 — 
tan _1 (— Q)/^] 1 ^ 2 , i.e. rms = R as Q — > oo). However, due 
to signal-to-noise considerations, NWD were unable to ob- 
tain an adequate fit above ~ 30 Hz. 

As the seven brightest observations of GX 339—4 pre- 
sented in NWD are so intrinsically similar, we decided to 
average these observations to form a composite PSD. The 
PSD, utilizing energy channels covering 0-21.9 keV (i.e., 
channels A-D from NWD summed) were calculated and 
averag ed together using the normalization of Leahy et al. 
( 1983j) (note that this normalization weights the signal PSD 
by the count rate of the observation, with the average count 
rates here ranging from 517-877 cps), the noise was sub- 
tracted (see Nowak et al. 1999a; and references therein), 
and then the PSD was renormalized to the normalization 



of Belloni & Hasinger (1990) using the mean count rate of 
the summed observations. For this normalization, integrat- 
ing over Fourier frequency yields the mean square variability 
relative to the square of the mean of the lightcurve . Erro r 
bars on the PSD were calculated as in Nowak et al. (1999a) 
We fit this PSD with a model that consisted of a zero 
frequency-centered Lorentzian (= A u /[1 + (///o) 2 ]), plus 
four additional QPO features of the form of eq. [j]. The two 
lowest frequency QPO features were constrained to be har- 
monics of one another (this allowed us to fit the "asym- 



© 0000 RAS, MNRAS 000, 000-000 



Are There Three Peaks in the Power Spectra of GX 339—4 an d Cyg X-l? 3 



POWER SPECTRAL DENSITIES FOR 
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Figure 2. Power spectral density (PSD) of GX 33 9-4 for the 
faintest observation presented in Nowak et al. (1999b). Thick solid 
line is the best fit model, whereas the dashed lines are the separate 
components of the model. The thin solid line is the residual noise 
level (see text) after subtracting the mean noise PSD. Lower panel 
shows the residuals from the fit. 



metric" nature of this feature that was commented upon 
in NWD), but otherwise all QPO amplitudes, widths, and 
frequencies were allowed to be fit parameters. Our best fit 
model is also presented in Fig. |lj The fit yielded \ 2 = 177 for 
75 degrees of freedom. The fit parameters (with subscripts 0, 
1, h, 2, 3, for the zero frequency-centered Lorentzian, lowest 
frequency QPO, its harmonic, and higher frequency QPOs) 
are presented in Table [l]. Error bars are 90% confidence level, 
i.e. A^ 2 = 2.71 for one interesting parameter. 

We have also applied the above fit to the faintest obser- 
vation (observation 5) from NWD. Here we are able to fit 
the zero frequency-centered Lorentzian, the low frequency 
QPO (plus harmonic), and the middle frequency QPO. The 
results are also presented in Table ^ and the data are shown 
in Fig. |2| The high-frequency QPO, if present, is essentially 
unconstrained (although in Table nlwe give the error bars for 
the "local minimum" in \ 2 )- Fig-0 a l so shows the residual 
noise level. This level corresponds to the expected ampli- 
tude of postive 1-er fluctuations above the mean value of 
the Poisson noise PSD, and therefore is indicative of the 
minimum PSD amplitude at which a signal can be detected 
(see Nowak et al. 1999a, and references therein). This resid- 
ual noise level, which scales inversely proportionally to the 
count rate and to the square root of the integration time, 
shows the difficulty of detecting the presence of a "third 
QPO feature" in the high frequency PSD of a single obser- 
vation. The composite PSD shown in Fig. [j] has a residual 
noise level « 9 times lower than that for the faintest obser- 
vation of GX 339-4. 



Nowak et al. (1999a) discussed observations of Cyg X-l in 
its low luminosity/spectrally hard state. In that work, we 
presented doublely broken power law fits to the PSD, with 
the low frequency PSD being essentially flat, the middle- 
frequency PSD being approximately oc /~ , and the high- 
frequency PSD being approximately oc f~ 2 . Although the 
fractional deviation of the data from the fits was quite small, 
the reduced \ 2 ranged from « 4-16 due to the excellent 
statistics achievable with RXTE. Here we refit the PSD with 
the functional form discussed in The fit parameters for 
two separate energy channels (0-4 keV and 14-45 keV; see 
Nowak et al. 1999a) are also presented in Table [j], and the 
fits are shown in Fig. |^. 

The Cyg X-l PSDs show subtle, but statistically sig- 
nificant, features. Broad peaks are seen at ~ 0.3 Hz, 2 Hz, 
10 Hz, and 50 Hz. The nearly flat portion of the PSD at 
Ss 1 Hz shows some sign of a peak. We use the same model 
that we fit to GX 339—4, namely a zero frequency-centered 
Lorentzian, a QPO and its harmonic, plus two additional 
QPO, and obtain reasonable fits. The reduced x 2 f° r these 
fits were significantly lower (A^ 2 > 150 for eight extra de- 
grees of freedom) than that for the double ly broken power 
law fits discussed in Nowak et al. (1999a). Note that the 
lowest frequency 'QPO' is identified with the portion of the 
PSD at £3 1 Hz, and its harmonic is identified as having both 
larger amplitude and greater Q value, in contrast to the fits 
to GX 339-4. 

Although these fits are a significant improvem ent ove r 
the doublely broken power law fits of Nowak et al. (1999a), 
there are large residuals, especially near 0.2 Hz. We thus 
considered another fit consisting of five QPO features plus 
one harmonic (or possibly sub-harmonic, as the QPO asso- 
ciated with the higher of the two related frequencies had a 
larger amplitude) . The results for these fits are presented in 
Table ^, and the fits are also shown in Fig. |j| This model 
resulted in significant improvements to the fits (A^ 2 ~ 50 — 
90, for four extra degrees of freedom); however, compared 
to the fits presented in Table |l| there was relatively little 
change in the parameters for the highest frequency QPO fit 
components. 

In these models, what was fit as a zero frequency- 
centered Lorentzian is now fit as a strong, broad QPO plus a 
weaker, narrow QPO. Whether these fit components should 
be regarded as separate physical features, or whether the 
improvement in the fit is merely an indication of the inade- 
quacy of the zero frequency-centered Lorentzian in describ- 
ing the low frequency PSD is, of course, unclear. Also along 
these lines, we note that although we have constrained the 
features at « 0.8 Hz and 1.6 Hz to be harmonics of each 
other, the fits are slightly improved (A^ 2 = 7.0 for the 0- 
4keV band and Ax 2 = 8.9 for the 14-45 keV band) if the 
frequencies are allowed to float freely. Note also that there 
are still positive residuals at 10 -3 Hz. The above points, 
as we discuss further below, serve to highlight the fact that 
when the PSD features are as broad and as subtle as they are 
in both Cyg X-l and GX 339—4 the intepretation of what 
is and is not to be considered a "QPO" becomes much more 
ambiguous. This needs to be borne in mind in the discussion 
below of the correlations between fit features. 
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Table 1. Parameters for PSD Fits. All frequencies arc units of Hz. 
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zfc-Lorcntz. 


QPOj 


QPO h 


QP0 2 


QPO3 


x 2 /DoF 
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A = 


„ ~~4-0 01 

u - z -0.01 
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_ . 4-0 01 

n 1 4^ 
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h = 
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/2 = 
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o.n+n? 


Ri = 


02+ 08 
u - uz -0.02 


Rh = 
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u - zu -0.05 
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Figure 3. Fourier frequency times the power spectral density (PSD) of Cyg X-l for the 14-45 keV band (left, middle) and the 0-4 keV 
band (right). Data are the same as presented in Nowak et al. (1999a). Thick solid lines are the best fit models from Table 1 (left) and 
Table 2 (middle, right), whereas the dashed lines are the separate components of the models. The thin solid lines are the effective noise 
levels after subtracting the mean noise PSDs. Lower panels show the residuals of the fits. 



Table 2. Parameters for PSD Fits to Cyg X-l Data. All frequencies are units of Hz. 
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x 2 /DoF 


R 
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QPOi 


29+ ' 02 
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QPO2 
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4 DISCUSSION 

The evidence for a third, high-frequency feature in the com- 
posite PSD of GX 339—4 is clear. Removing this feature, 
the x 2 increases to 442. Here we should note that to some 
extent the fits are affected by the fact that separate ob- 
servations with slightly varying PSD properties have been 
averaged together. Some broadening of the "QPO features" 
is expected as a result, and the Q values for features fit 
to the summed PSD should be considered as lower limits. 
Likewise, the fitted frequencies and amplitudes of the vari- 
ability features should be viewed as indicative of an average 
value, but not strictly applicable to any of the individual 
PSD used in the combined observation. Furthermore, their 
error bars should be viewed as lower limits. Along these 
lines, we note that although the presence of a third feature 
is indicated, it is almost as well-modelled (x 2 = 207) by a 
zero frequency-centered Lorentzian (oc [1 + (///&) 2 ] -1 ) with 
a break frequency of ft = 23 ± 3 Hz. 

The question arises as to what extent this third variabil- 
ity feature might be an artifact of averaging together seven 
separate observations. This is of some concern; however, we 
note that in terms of well-measured properties, the seven 
individual PSDs are very similar to ea ch other. Using the 
normalization of Belloni & Hasinger (Belloni & Hasinger 



199C), at frequencies < 10 Hz (where all seven individual 



PSDs have good statistics) the variance of the noise sub- 
tracted PSDs is < 25%. Specifically, 



1/2 



Ps(fT 



[Ps(f) - P,(f)] 2 



< 0.25 



(2) 



at all frequencies / < 10 Hz, where Pi(f) are the noise sub- 
tacted PSD of the individual observations, and P s {f) is the 
noise subtracted PSD of the summed observation. The vari- 
ance of the mean is lower by a factor of \/6 and therefore 
is ^ 10%. This result is essentially unchanged if we weight 
the observations by count rate. The fact that some varia- 
tion exists from observation to observation means that the 
error bars that we used should be viewed as lower limits, 
and the derived \ 2 should be viewed as upper limits. Also 
as regards possible systematic variations, NWD showed the 
broad features at ~ 0.35 and 2.5 Hz have fitted frequencies 
that ranged only from ~ 0.3-0.4 Hz and 2-3 Hz, respectively, 
among the individual observations and among separate en- 
ergy bands within each of these observations. Given the sim- 
ilarities among the PSDs, it would be unusual for the strong 
third variability feature present in the combined observation 
to be an artifact of the PSD averaging. (We note that the 
third variability feature discussed below for Cyg X-l was for 
a single, short observation.) 

Taking all three features in GX 339—4 as real and not an 
artifact of the PSD summation, these features are relatively 
broad (even considering possible systematic broadening due 
to averaging) compared to high frequency QPO seen, for ex- 
ample, in neutron stars. Statistically, they all have Q > 
at a high significance level; however, they are all far from 
narrow, periodic features. The QPO at the higher frequency 
end of the correlation presented by PBK all tend to be much 
narrower features (in terms of FWHM) than those presented 
here. Nevertheless, if we take these features and identify the 
low-frequency QPO with the HBO, the middle-frequency 



QPO with the lower-frequency kHz QPO, and the high- 
frequency QPO with the upper-frequency kHz QPO, then 
they approximately fall along the correlations suggested by 
PBK. We show the above data overlaid on the suggested 
correlation of PBK in Fig. This is true for the composite 
PSD as well for the PSD of the faintest observation from 
NWD, keeping in mind that for this latter observation the 
error bars on the highest frequency feature are for the lo- 
cal minimum and that the presence of this latter feature is 
not strongly constrained. We also note that the correlation 
between the frequency, fo, of the zero frequency-centered 
Lorentzian and the frequency, /i, of the lowest frequency 
QPO is in the same sense as the break frequency /QPO fre- 
quency correlation discussed by Wijnands & van der Klis 
( jjgj )- 

The fact that the "HBO" and "lower-frequency kHz 
QPO" of GX 339—4 appear to lie approximately along the 
suggested correlation was already noted by PBK. What is 
new here is that there appears to be a third QPO that ex- 
tends the lower- frequency kHz QPO/upper- frequency kHz 
QPO correlation downward by nearly two orders of mag- 
nitude in Fourier frequency. The high-frequency QPO of 
GX 339—4 appears to lie slightly below the extrapolation of 
the lower- frequency kHz QPO/upper- frequency kHz QPO 
correlation; however, assuming the correlation to be oc f a , 
a slope change of 8a ~ —0.2 is all that is required between 
the low-frequency GX 339—4 point and the high-frequency 
trend. Furthermore, some theories (e.g., the relativistic pre- 
cession theory of Stella & Vietri 1998; Merloni et al. 1999; 
Stella et al. 1999) suggest that one should not expect a single 
power law over the whole range of the putative correlation. 
Specific predicted frequency correlations for th e rela tivistic 
precession theory can be found in Stella et al. (1999). 



The features fit to the Cyg X-l PSD are slightly more 
problematic to fit onto this trend. In light of the fact that 
the best fit consists of five QPO features plus one harmonic, 
there are the questions of which fit components are to be 
considered as "true QPO" (as opposed to merel a fit arti- 
fact) and which are to be associated with which feature in 
the putative correlation. In Fig. ^j, we have taken the fits 
of Table ^| and associated QPO5 with the upper-frequency 
kHz QPO, QPO4 with the lower-frequency kHz QPO, and 
QPO3 — the sub- harmonic to the stronger 0.16 Hz feature — 
with the HBO. With these associations, the features fit to 
the high energy (14-45 keV) PSD of Cyg X-l agree with 
suggested trend, whereas the features fit to the low energy 
(0-4 keV) PSD of Cyg X-l deviate from the trend, mostly 
due to the extremely low frequency of the "lower-frequency 
kHz QPO" . We note, however, that the frequency-break of 
this component occurs at roughly the same location in the 
low and high energy PSDs. The fact that the feature in the 
low energy PSD has an extremely low Q ~ 0.2 leads to a 
very low fit-frequency in order to yield a break in the same 
location as the somewhat narrower feature (Q « 0.6) in the 
high energy PSD. If one performs a joint fit to the 0-4 keV 
and 14-45 keV PSDs, constraining the QPO frequencies (but 
not the amplitudes or widths) to be the same in both bands, 
a reasonable fit is obtained with x 2 = 156 for 111 degrees 
of freedom. The resulting fit frequencies are comparable to 
the frequencies fit to the 14-45 keV band in Table ^. 

It is difficult to test the correlations suggested by PBK 
since most of the Q values for the fitted features are very 
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Figure 4. The QPO frequencies presented here (filled circles are 
for GX 339-4, filled squares are for Cyg X-l), with 90% CL error 
bars, overlaid on the sa me graph as presented in Psaltis, Belloni, 
&c van der Klis (1999). Open diamonds correspond to neutron 
star candidates, filled diamonds to black hole candidates (error 
bars are not shown for either). Note: QPO identified as "normal 
branch oscillations" are also shown on this plot (see PBK). 



low, making the features very subtle. Another, more serious 
problem, is the question of identification of the features. One 
alternative fit to the Cyg X-l PSDs would be to place a zero 
frequency-centered Lorentzian with break frequency ~ 0.01 
to remove the low frequency residuals seen in Fig. |§|, and 
then to identify QPOi and QPO2 as the HBO and a har- 
monic, QPO3 and QPOh as the lower-frequency kHz QPO 
and a harmonic, QPO4 as the upper- frequency kHz QPO, 
and QPO5 as a new, "fourth QPO". [The theoretical ex- 
istence of a fourth such frequency has been suggested by 
Psaltis & Norman (2000), for example.] These frequencies 
would also fit on the suggested trend of PBK, except now 
QPO4 in the low energy PBK would be too low a frequency 
to fall on the trend for the upper- frequency kHz QPO. At 
present, there is no unambiguous, rigorous method of asso- 
ciating a fitted feature with the suggested trend of PBK. 

Even given these caveats regarding correlations between 
the various fit-components, it is tempting to associate each 
fit-component of the PSD with a "resonance" in the variabil- 
ity properties of the accretion flow (see Psaltis & N orman 
2000). This is counter to the m odels of Kazanas et al. (1997 ) 



and Poutanen & Fabian (1999), for example, that essentially 
postulate a single "response" for the variability properties 
of the accreting system. NWD phenomenologically elabo- 
rated upon the concept of each PSD fit-component corre- 
sponding to a separate physical mechanism by further pos- 
tulating that each of these system "responses" had a sep- 
arate "driver", uncorrelated with the variability-drivers of 
the other PSD components^]. The net Fourier frequency- 

* Even if a system has a set of independent responses to a source 
of input fluctuation or noise, for example the accretion disk re- 



dependent phase lags (or, equivalently, time lags=phase 
lags/27t/) and coherence function [j 2 {f ), a measure of the 
degree of linear correlation between two lightcurves; see 
Bendat & Piersol 1986, Vaughan & Nowak 1997] between 
soft and hard variability were then given by a combination 
of the phase lags and coherence functions for each individual 
fit-component of the PSD (see NWD, §4 and Fig. 8; Vaughan 
& Nowak 1997). 

NWD noted that for the case of GX 339-4 wherever 
one PSD fit component dominated (for example, near 0.3 Hz 
in Fig. |lj where the low-frequency QPO dominates) , the co- 
herence function would be near unity and there would be an 
approximately flat Fourier phase lag shelf when comparing 
soft and hard variability. Wherever two PSD fit components 
crossed one another, there would be a slight dip in the co- 
herence function and a transition from one characteristic 
Fourier phase lag to another (see Fig. |). Fig. g shows that 
at any given Fourier frequency, typically two PSD fit com- 
ponents dominate. Ignoring the other two fit components, at 
any given frequency we can then calculate the phase lag be- 
tween soft and hard X-ray variability for the remaining two 
PSD fit components by simultaneously fitting the measured 
net phase lag and the measured net coherence function^. 

We have carried out such a procedure for the sec- 
ond brightest GX 339—4 observation (which showed the 
strongest QPO at « 0.3 Hz) presented in NWD. We have fit 
the 0-21.9 keV PSD of this observation with the same type of 
model as discussed in Table For this non-composite PSD 
fit, constraints are again weakest for the "third QPO" fit- 
component. Furthermore, we have assumed this PSD shape 
and amplitude for both the 0-3.9 keV and 10.8-21.9 keV 
PSDs. (Fitting each energy band individually, no constraints 
can be made on the "third QPO".) We then followed the 
procedures outlined in §4 of NWD (see specifically eqs. 4 
and 5) for fitting the phase lag and coherence function, ex- 
cept that instead of fitting four constant (as a function of 
Fourier frequency) phases over the entire data set, we fit two 
phases at any given Fourier frequency. Specifically, at any 
given Fourier frequency, we fit the phase lags for the two 
strongest PSD fit- components, assuming the other two PSD 
fit components to be negligible. Error bars do not account 
for uncertainties in the PSD fit parameters, and only repre- 
sent that portion of the error due to the uncertainty in the 
measured phase lag and measured coherence function. The 



sponses discussed by Psaltis & Norman ( t200(l| ) . the net outputs 
will be perfectly correlated (i.e., have a coherence function of 
unity) if they are responding to the same source of noise fluc- 
tuations "This point is discussed i n furt her detail by Bendat & 
Piersol ( 1986[ ), Vaughan & Nowak ( 1997 ), and NWD. 
t Note that there is still an ambiguity in the values of the 
Fourier phases, as there are two independent solutions for a two- 
component fit to the measured phase lag and coherence function. 
Specifically, see Fig. 8 of NWD, where reflecting the phases of 
the individual PSD components through the net measured phase 
lag yields and equally valid solution. In our solution that follows, 
at low Fourier frequency we choose the phase lag for the zero 
frequency-centered Lorentzian to be the value closest to the net 
measured phase lag. Then, progressing to higher Fourier frequen- 
cies, we choose solutions that are most nearly continuous with 
the lower frequency values. 



© 0000 RAS, MNRAS 000, 000-000 



Are There Three Peaks in the Power Spectra of GX 339—4 an d Cyg X-l? 7 



.10 ^ 

G .08 -- 

V 06 ■- 

± .04 r 

.02 I- 



1.0 



| 0i- o O 

M <■> 

i—l 

<D 
aa 
P3 



*<[> 

O 

<J> ♦ ♦ * 



□ - 



o.o -^-hr^ 



-0.5 



10-' 



1. 

Frequency (Hz) 



10. 



Figure 5. Bottom Panel: A possible decomposition of the net 
Fourier frequency-dependent phase lags between soft (0— 3.9keV) 
and hard (10.8—21.9 keV) X-ray variability (filled circles) into indi- 
vidual components associated with the PSD fit components shown 
in Fig. |lj. Triangles correspond to the zero frequency-centered 
Lorentzian, diamonds to the low-frequency QPO (including har- 
monic), clear squares to the middle-frequency QPO, and filled 
squares to the high-frequency QPO. Positive values correspond 
to the hard X-ray variability lagging the soft. Top Panel: The 
coherence function, 7 2 (/), associated with the net phase lag. 



results for this decomposition (which, as we discuss above, 
is not unique) are presented in Fig. |E| 

Fig. ^ is presented more for illustrative purposes rather 
than as a completely rigorous and definitive fit. The things 
to note about this possible decomposition are the follow- 
ing. If the QPO fit components each represent a sepa- 
rate, independent physical process, then the deviations from 
unity coherence can be understood within the context of 
such a model as well. At low frequency the net phase/time 
lags are most characteristic of the zero frequency-centered 
Lorentzian, at middle frequencies they are most represen- 
tative of the low frequency QPO, and at moderately high 
frequencies they are most representative of the middle fre- 
quency QPO. At w 0.3 Hz, there is a recovery of the co- 
herence function to near unity values that would be asso- 
ciated with the low frequency QPO becoming dominant in 
the PSD. Although it is difficult to accurately measure the 
high frequency phase/time lags and coherence function, the 
fact that the middle and high frequency fit components to 
the PSD have comparable amplitudes at / <; 10 Hz could ex- 
plain the strong loss of coherence at these frequencies (which 
is characteristic of Cyg X-l as well; see Nowak et al. 1999a). 
Although the net phase lag indicates that hard photons lag 
the soft photons, it is possible for portions of the individual 
variability components to show exactly the opposite behav- 



5 SUMMARY 

The main results of this work can be summarized as follows. 



• A composite power spectrum of GX 339—4 is well- 
fit by a model that consists of a zero frequency-centered 
Lorentzian, a moderate width (Q ~ 1-2) quasi-periodic 
oscillation and its harmonic, plus two additional, broad 
(Q ~ 0.5-1) QPOs. 

• The frequencies of these three QPO components ap- 
proximately agree with the correlations among horizon- 
tal branch oscillations/lower- frequency kHz QPOs/upper- 
frequency kHz QPOs suggested by Psaltis, Belloni, & van 
der Klis fll999| ). 

• Fits to the PSD of Cyg X-l are more problematic. Mul- 
tiple QPO components provide the best fits to the PSD as 
a function of observed energy band. Depending upon the 
identifications made, the high energy band PSD fits can be 
made consistent with the corre lations suggested by Psaltis, 
Belloni, & van der Klis fll999| ). The low energy band PSD 
fits show deviations from these correlations, although a joint 
fit of the low- and high-energy PSD does allow for a single 
set of frequencies consistent with the suggested correlations. 

• One can construct a phenomenological model where the 
net Fourier frequency-dependent phase lag and coherence 
between soft and hard X-ray variability can be expressed 
as the combination of phase lags and coherences from the 
individual fit components to the PSD. 

• If such a model for the phase lags and coherences is 
correct, then although the net phase lag indicates that soft 
variability leads the hard variability, the trend may be ex- 
actly opposite for some of the individual variability compo- 
nents. A requirement of this model is that although the in- 
dividual components have coherent variability between soft 
and hard photons, the variability components are incoherent 
with one another. This not only requires that each compo- 
nent be representative of a different "response" within the 
accretion flow, but also that each component has a separate 
variability "driver" 

The correlations between the different fit-components 
seems at least to be very suggestive. It is tempting to as- 
sociate all of these quasi-periodic features observed in both 
neutron star and black hole systems with a common set 
of phenomena intrinsic to the accretion flow itself, and not 
to the properties of the compact object such as the pres- 
ence of a surface or a magnetic field. Such may not be 
the case in reality, however. The presence of what appear 
to be four distinct features (the zero frequency-centered 
Lorentzian and three broad, quasi-periodic features) in the 
PSD of GX 339—4, coupled with suggestive evidence from 
the Fourier frequency-dependent phase lags and coherence 
function (see also NWD), at the least seems to be arguing 
for models with multiple, distinct sources of variability (e.g., 
Psaltis & Norman 2000), and against models with a single 
"type" of variability (e.g., Kazanas et al. 1997, Poutanen & 
Fabian 1999). 
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